Theoretical Study of Pressure Effect on TDAE-C 60 
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We have theoretically studied pressure effects on molecular ferromagnet Ceo complexes with tetrakis 
(dimethylamino) ethylene (TDAE), particularly the pressure- induced depression of the Curie tem- 
perature. The observed behavior is well simulated with our model based on a charge transfer induced 
intramolecular Jahn- Teller distortion and an intermolecular cooperative Jahn- Teller interaction. We 
emphasize that the theoretical simulation is carried out with reasonable parameters known for C^ 
complexes. It is concluded that the enhancement of crystal field at Ceo site due to increasing pressure 
causes the depression of Curie temperature. 



In 1991, Allemand et al. reported the ferromagnetic 
behavior of tetrakis (dimethylamino) ethylene (TDAE) 
-Ceo with Tc—16 K jy. It has been attracting many 
scientists due to its having the highest Curie temper- 
ature among pure organic molecular ferromagnets and 
its unusual magnetism. Allemand et al. proposed 'soft 
ferromagnetism' with no hysteresis in the M-H curve jj] . 
Tanaka et al. Q] and Blinc et al. || have observed behav- 
iors like a 'superparamagnetism' of spin clusters consist- 
ing of hundreds of spins in the magnetization and the pro- 
ton nuclear magnetic resonance (NMR) measurements 
below To, respectively. Venturini et al. suggested 'spin 
glass model' from analysis of the electron spin resonance 
(ESR) lineshape Q|. In 1997, we proposed a model for 
clarifying the properties of this system || . This model is 
introduced in order to clarify the origin of intermolecular 
ferromagnetic coupling between Cgo's, and it is based on 
the orbital ordering of unpaired electrons on Ceo's due to 
the adjacent alignment of the Jahn- Teller distorted Cgo's. 
The spin-glass-like behavior and superparamagnetism of 
the spin clusters may also be explained qualitatively in 
this framework ||. 

Very recently, Mizoguchi et al. observed the pressure 
dependence of the Curie temperature (Tc) in TDAE-Ceo, 
as shown in Fig. |l| . It appears that Tq is parabolically 
depressed upon the application of pressure. The purpose 
of this paper is to provide a quantitative understanding of 
this behavior in the framework of the model mentioned 
above. It should be noted that the observed value is 
simulated quantitatively with the reasonable parameters 
known for the Gqq molecule. 

First, our orbital ordering model is briefly introduced. 
With respect to this model, the magnetic interaction be- 
tween local spins on Ceo's are discussed. In TDAE-Cgo, 
both TDAE and Cgo are regarded to have an unpaired 
electron since one electron transfers from TDAE to Ceo- 
However, it has not yet been clarified whether or not bare 
spin moments exist on TDAE molecules. On the other 
hand, unpaired electrons on Cgo molecules are considered 
to play a crucial role in bulk ferromagnetism. Figure ^ 
shows an example of a molecular arrangement which is 
likely to cause three-dimensional ferromagnetic order. 



The lowest unoccupied orbitals (LUMO) of Ceo in Ih 
are triply degenerated with t\ u symmetry. When degen- 
erated orbitals are partially occupied, the molecule is dis- 
torted in order to stabilize one orbital, which is known 
as the Jahn- Teller (JT) effect. We assume the T>2h struc- 
ture for the C^q anion in this study. We performed a geo- 
metrical optimization of single a C^~ anion in a previous 
study. The optimized structure resembles a rugby ball 
whose elongated axis lies along one of the three symme- 
try axes in the T>2h structure. The atomic displacement 
from the ideal icosahedral is small, i.e., 0.01 A at most. 
The three t lu states, LUMO^, LUMO y , and LUMO z , 
which are degenerate LUMO's of Ceo with 1^ symmetry 
before the charge-transfer, each split into three states. If 
the elongated axis is the x-axis, LUMOz has the lowest 
energy, and the other two orbitals, LUMO y and LUM0 2 , 
with higher energies are almost degenerate. It should be 
noted that the charge density is not spherical but takes 
large values along a belt surrounding the elongated axis 
(see Fig. |). 

Magnetic interactions between distorted molecules de- 
pend on their alignment. We concluded that the ferro- 
magnetism is realized in the structure illustrated in Fig. ^ 
due to periodic molecular distortion In a C^ ar- 
ray along the c-axis in this structure, the elongated axes 
of two nearest C^~ 's are perpendicular to each other. 
In this case, the orbitals of unpaired electrons align in 
an alternating manner, for example, LUMO x , LUMO y , 
LUMO x ,... Such a system is called an orbital ordering 
system and an interaction acting to favor an alternating 
alignment is called a cooperative JT interaction. It is 
known that the ferromagnetic coupling between neigh- 
bors is preferred in such a orbital ordering system. It 
should be noted that the intermolecular transfer occurs 
only between the same kind of orbitals if the crystal field 
at the C^g's has perfectly orthorhombic symmetry. It is 
a key point in favoring the intermolecular ferromagnetic 
interaction. 

It is expected that the magnetic interactions between 
Ceo's lying along the c-axis are the strongest among the 
intermolecular magnetic interactions, because their dis- 
tance is shorter than that in the ab-plane by 0.3 A. Elec- 
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trical transport along the c-axis is also observed about 
10 times larger than that along the a-axis J^]. In this pa- 
per, we assume that the system is quasi-one-dimensional 
along the c-axis. The extended Hubbard Hamiltonian for 
the one-dimensional chain with LUMO x 's and LUMO y 's 
as the basis functions is shown as follows: 
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where cj^^c^o-) represents the creation (annihilation) 
operator of /z-orbital with cr-spin at the i-th Cgo molecule 
and rii^a is its number operator. The lines represent 
the orbital energy, the intermolecular transfer energy 
between different orbitals, that between the same or- 
bitals, the intra-orbital Coulomb energy, the inter-orbital 
Coulomb energy and the S Z S Z part of the exchange in- 
teraction, and the S+S- part of exchange interaction, re- 
spectively. The intra-orbital and inter-orbital Coulomb 
energies are assumed to be the same. LUMO z 's are not 
included because they are not occupied on any Cgo's in 
the molecular alignment shown in Fig ||. It should be 
noted that a small transfer energy exists even between 
different orbitals related to the third term because of a 
small deviation from the perfect orthorhombic symmetry 
of the crystal. 

We suggest that the unpaired electron can be regarded 
to exist locally on each molecule since the Coulomb en- 
ergy is sufficiently larger than the transfer energy. The 
energies ti and U are estimated to be about 0.05 eV and 
0.6 eV, respectively P,10[. Experimentally, the Curie- 
Weiss like behavior is also observed in the magnetic sus- 
ceptibility even under the application of pressure |?J . For 
a localized system, the Hamiltonian ([!) can be trans- 
formed into an extended Heisenberg Hamiltonian. 
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where J\ and J2 represent intrachain exchange inter- 
actions between nearest neighbors and second nearest 



neighbors, respectively. The magnitude is calculated 
from a second order perturbation of the transfer energy 
and the maximum terms of a fourth order perturbation. 
Details of the transformation are reported for another 
orbital ordering system, K2CUF4 [jy]. From this Hamil- 
tonian, the Curie temperature Tq is derived with mean 
field theory as 
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where we assume the existence of interchain ferromag- 
netic couplings J3, and k represents the Boltzmann fac- 
tor. £(=1/2) and J(0)(= 2Ji + 2J 2 + 4J 3 ) indicate the 
value of a single spin and the q=0 component of the 
Fourier transformation of exchange interactions J(q), re- 
spectively. It should be noted that the mean field theory 
generally overestimates the Curie temperature Tq. Here 
we derive the pressure dependence of TDAE-C60 based 
on Hamiltonian |^. It is assumed that the transfer energy 
linearly depends on pressure as t a = t° a +pt' a (a = £,g). To 
describe the difference of the Curie temperature derived 
by mean field theory T^f F from that observed experi- 
mentally Tq* , a reduction parameter a is introduced as 
T<? x = aT™ F . 

The solid line shown in Fig. [j] represents the simulation 
result of the pressure dependence of Tq with £°=0.065eV, 
i°=0.0035eV, 4=0.001eV/kbax, ^=0.00243eV/kbar, 
C/=0.55eV, J=0.09eV, A=0.15eV, J 3 =3K, and a=0.75. 
It is found that the theoretical result simulates the ob- 
servation well. We emphasize that the parameters used 
in the simulation are reasonable. The intramolecular pa- 
rameters U, J, and A are almost the same as the values 
estimated by Suzuki and Nakao [|l0| . The report of ab- 
initio calculations for fcc-Ceo by Saito and Oshiyama is 
used for the determination of the intermolecular transfer 
t° §]. The width of the LUMO band is calculated to be 
about 0.5 eV in fcc-Ceo and the intermolecular transfer 
energy can be roughly estimated as 0.04eV. It is in good 
agreement with our parameter t®. We have no quantita- 
tive information of the intermolecular interaction J3 . It 
is considered that 3K may be appropriate for the value of 
J3 because the interchain distance is also longer than the 
intrachain molecular distance by about 0.3 A. In conclu- 
sion, the parameters used in our simulation is found to be 
quite reasonable for TDAE-Cgo ■ Tanaka et al. calculated 
the dependence of the intermolecular magnetic coupling 
between Cgo's on their orientation with a semi-empirical 
approach without the JT distortion jf^] . They reported 
that even the strongest ferromagnetic coupling is very 
small, about 0.03 K. Therefore, it is very important that 
the high Tq of this material can be explained by the JT 
distortion of Cgo's with the reasonable parameters. 

Figure ^ shows the pressure dependence of the inter- 
molecular exchange interactions Ji and Ji as well as their 
sum. It is found that the negative parabolic shape in 
the pressure dependence of Tq is almost determined by 
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Ji although its magnitude is weakened by the antiferro- 
magnetic interaction J 2 . The negative dependence of Tq 
on the application of pressure is formed by a coefficient 
of the parabolic term in the pressure dependence of J\ 
shown in eq.^. With the parameters used in the simu- 
lation, it is expected that Tq is decreased by applying 
pressure in the case of t\ < 12ig. The parameters for the 
simulation satisfy this condition. The crystal structure 
of TDAE-Ceo is slightly different from the orthorhombic 
one although the Ceo molecule has complete orthorhom- 
bic symmetry. Therefore, the enhancement of the crystal 
field at Ceo's site due to the application of pressure makes 
t g significantly larger. In fact, it is suggested from an ab 
initio calculation that t' g and t' e have the same order in 
another orbital ordering system, K2CUF4, in which local 
crystal field at the Cu site also has a lower symmetry 
than that of JT distorted molecules [jll) . 

A helical magnetism can be realized if the second 
neighbor coupling J 2 is sufficiently large 



amine the stability of ferromagnetism against the heli- 
cal magnetism. According to the mean field theory, the 
magnetic order with the wave number Q is realized if the 
Q-component is the largest in the Fourier form of an ex- 
change interaction J(q). In the case of J1/IJ2 > 4, a 
ferromagnetic ordering is the most stable. Otherwise a 
helical magnetism is realized. It is evident in Fig. || that 
helical ordering is expected above 6 kbar. It is not clar- 
ified at the present time whether the magnetic ordering 
at p=7A kbar is ferromagnetic or helical ordering with a 
long wavelength. 

Let us discuss the dimensionality of TDAE-Cgo and the 
origin of the interchain ferromagnetic interaction. We as- 
sume the one-dimensionality of this system in our model. 
It is also indicated experimentally with the measure- 
ment of the electrical conductivity ||. Although Blinc 
et al. reported that TDAE-Cgo is the isotropic ferro- 
magnet in the temperature dependence of the spin-wave 
resonance in ESR fl3|j , the temperature dependence of 
the magnetization of the quasi-one dimensional system 
with J1/J3 ~ 6 is not significantly different from that 
of an isotropic system. Therefore, the assumption of 
quasi-one-dimensionality does not conflict with the ob- 
servations by Blinc et al.. The origin of interchain fer- 
romagnetic coupling can be explained with the orbital 
ordering structure shown in Fig. |[ In this structure, 
the ferromagnetic ordering is favored not only inside the 
chain but also interchain coupling || . The quasi-one di- 
mensionality can be explained with this mechanism since 
the magnitude of magnetic interaction depends on in- 
termolecular transfer energy. The interchain distance is 
about 0.3 A longer than the intrachain molecular dis- 
tance. Another model is the superexchange mechanism 
by way of TDAE molecules. The status of unpaired elec- 
trons on TDAE+ cations is an open question but it is 
likely that their molecular spins form pairs The 
reason for electron pairing on TDAE should be clarified 
in order to discuss the interchain magnetic interaction. 

In this paper, we assumed that one two-fold axis is par- 



allel to the c-axis. However, the constraint for realizing 
the ferromagnetic interaction along the c-axis is looser 
in reality. The suppression of the electronic transfer in- 
teraction between singly occupied molecular orbitals of 
neighboring Cgo anions is the most important. Conse- 
quently, our model requires that the elongated axes of 
neighboring Ceo's be oriented perpendicular or parallel 
to the c-axis and also perpendicular to each other. A 
structure proposed with X-ray analysis jl5| is feasible 
with this constraint. 

Recently, additional collateral evidence has been re- 
ported. Kambe et al. also observed a structural phase 
transition at 180 K with X-ray diffraction measurements 
and they have reported the possibility of cooperative JT 
ordering in the low-temperature phase fll?! . Clear exper- 
imental evidence of D2h JT distortion in monoanion-Cgo 
was also reported on a single crystal of a model com- 
pound of monoanion-Cgo, [As(C6H 5 )4] 2 C6oCl |17|. Fur- 



1 1 1 - We ex- thermore, in another C^~ compound, an X-ray analysis 
also suggests the stabilization of the static JT distortion 
by symmetry lowering of the crystal caused by rotational 
ordering Jl8[ . For complete confirmation of our model, 
we urge to perform the neutron scattering observation. 
In neutron scattering measurements, we can obtain in- 
formation on the spin density. The spin distribution on 
C^q's is not spherical but takes large values along the belt 
around the elongated axis, as shown in Fig. ^. It has a 
much larger spatial difference than the lattice distortion. 
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FIG. 1. The pressure dependence of Curie temperature in 
TDAE-Ceo- Closed circles represent observed value. Solid line 
indicates the simulation result based on the orbital ordering 
model. 




FIG. 2. An example of possible JT-distorted crystal struc- 
tures which are likely to cause a three-dimensional ferromag- 
netic ordering. Elongated spheres indicate JT-distorted Ceo's 
and thick lines express one-dimensional chains. Note that dis- 
tortion is exaggerated. In this structure the elongated axes 
of Ceo's are perpendicular to each other for the interchain 
nearest neighbor Ceo's as well as for the intrachain nearest 
neighbor Ceo's. The gray belts around Ceo's schematically 
represent the distribution of a unpaired electron. 
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FIG. 3. The pressure dependence of exchange interac- 
tion between intrachain nearest neighbors Ji , between second 
nearest neighbors J2, and their sum. 
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